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Abstract. We provide a proof to a recent conjecture by Forrester (2014, J. Phys. 
A: Math. Theor. 47 , 065202) regarding the algebraic and arithmetic structure of 
Meijer G-functions which appear in the expression for probability of all eigenvalues 
real for product of two real Gaussian matrices. In the process we come across several 
interesting identities involving Meijer G-functions. 
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1. Introduction 

Product of random matrices and the associated eigenvalue spectra exhibit a number of 
interesting properties and hnd concrete applications in varied helds of knowledge. Their 
study goes back to as early as 1950’s where the focus was on exploring the behavior 
of dynamical systems, and the accompanying questions related to stochastic differential 
equations and Lyapunov exponents pg. In the last few years there has been a revival in 
interest in their investigation because of their fascinating integrability properties [51-fT^ 
and identihcation of new problems where they can be applied, such as random graph 
states 1201. combinatorics izg, quantum entanglement ra and multilayered multiple 
channel telecommunication |9l|22]. 

In a recent study by Lakshminarayan m it was shown that the question concerning 
optimal quantum entanglement is intimately related to the probability that the product 
of two real 2-dimensional Gaussian random matrices (real Ginibre matrices) has real 
eigenvalues. This problem was exactly solved in [H] and extensive numerical exploration 
was conducted for the product of arbitrary number of matrices of higher dimensionalities. 
A complete analytical solution was provided by Forrester in [15] where he found explicit 
determinantal expressions for the probability of all eigenvalues real for product of any 
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number of Gaussian matrices of arbitrary dimension. Several related asymptotic results 
were also derived. In the finite dimensionality case the determinants were found to 
involve certain Meijer G-functions. Meijer G-functions have also appeared as correlation- 
kernels in product of Ginibre matrices or of truncated unitary matrices pf flOlfTSlIT^ . 
Bures and Gauchy two-matrix models Hi, and very recently in the results for product 
of a Wigner matrix and a Wishart matrix |23]. We give a brief introduction to Meijer 
G-functions in the next section. 

One of the intriguing observations made in ca involved the Meijer G-function 
5/2-j, 5/2-j, 2 
1, 1 -j- /l, \ k 

of two Gaussian matrices. Based on high-precision numerical computation it was 
conjectured that these are rational multiples of for positive integers j, k. We provide 
here a proof to this conjecture by deriving an exact and simple expression which clearly 
demonstrates why this is the case. In the process we also come across some interesting 
identities involving Meijer G-functions. 


^3,2 

'-^ 3,3 


which appeared as the kernel in the case of product 


2. Exact evaluation of a class of Meijer G-functions 


Meijer G-function is defined in terms of Mellin-Bernes type contour integral 

ds IIT=1 - s) Il]=i r(l - aj + s) 
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ni=m+i r(i - bj + s) Yrj=n+i r(ai - s)' 

where an empty product is interpreted as 1. The integer indices satisfy 0 < m < 
q,0 < n < p, and the parameters a’s, 6’s may be real or complex such that no pole of 
UT=inb,-s) coincides with any of the poles of 11^=1 r(l — Oj -|- s). Also, the contour 
C is such that it separates the poles of — s) from those of nj=i ^(1 — cij + s); 

see [211123] for details. One of the fascinating properties of the family of Meijer G- 
functions is its closure property under differentiation as well as indefinite integration. 
Moreover, majority of the established elementary functions and special functions can be 
represented in terms of Meijer G-functions. A few examples are 
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Uz) = G]fl I ^ 


Z//2, —V12 


TT / \ ^ 

- < aigiz) < 


( 6 ) 


Kiz) = I 


1/2, 1/2 

0 , 0 


(7) 


The special functions on left hand side of Eqs. ([5])-([7]) are upper incomplete Gamma 
function, Bessel function of the hrst kind, and complete elliptic integral of the hrst 
kind, respectively. 

For our purpose we begin with the identity 


■y3,2 


Gfs I ^ 


( 8 ) 


Ui, 02 , c \ _ 2,2 ( ®2 

The above result follows because of the presence of common parameter c in the upper 
and lower sets of parameters in the Meijer G-function on the left-hand side pTfEB] . The 
Meijer G-function in the above equation has a representation in terms of regularized 
Gauss hypergeometric function 


G: 
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— r(i — oi -|- 6i)r(i — 02 -i- 6i)r(i — oi -i- 62)7(1 — 02 -1- 62) 


2 Ti (1 — tti bi, 1 — 02 + 61 ; 2 — oi — 02 + 61 + 62 ; 1 — z). (9) 

We note that the regularized Gauss hypergeometric function is related to the 
usual Gauss hypergeometric function by the simple relation 2 Fi{ai,a 2 ] /3] z) = 
2 Fi{ai,a 2 ; /3] z)/T{/3). We now use one of the derivative identities for Meijer G- 
function 
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which leads, when applied recursively, to the following expression: 
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This is an interesting result which can be directly used along with (|9]) to demonstrate the 
truth of the conjecture in na. However, we make further progress and insert relation ([9]) 
in the above equation and use the following results to carry out the n-fold differentiation: 
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■2Fi{a,/3;^;l-z) = {-l)"^{a)m{/3)m2Fi{a+m,/3+m;^+m;l-z).{14:) 


Here 


^ ^ j represents the binomial coefficient and {a)n = r(a + n)/r(a) is 

the Pochhammer symbol. Incorporating all these results and carrying out some 
simplihcation, we arrive at the following identity: 


Gg I ^ 


Oi, 02 , C 


= r(l - Oi + 62)r(l - 02 + &2) 
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Substituting z = 1 and using 


2 Fi(a,/ 3 ;r, 0) = 


r(7)’ 


we obtain 
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Equations flT^ and ffT7|l are two key contributions of this work. We use ffT7|l in the next 
section to prove the result predicted in [15] . 


3. Probability of all eigenvalues real for product of two real Gaussian 
matrices 

Let us consider two iV x Wdimensional real Gaussian matrices X and Y which have 
the associated joint probability measure 

Ar2 


r{X,Y)d[X]d[Y] = 


27r 


e-h^dxx^+YYY^^X]d[Y]. 


(18) 


Here tr represents trace, and T stands for transpose. Also d[X] refers to the product of 
differentials of all matrix elements of X. Similar dehnition is to be understood for d[Y]. 
Forrester has derived an exact determinantal expression for the probability of all N 
eigenvalues real for product of m number of x Wdimensional standard real Gaussian 
matrices na. We quote here the result for m = 2, i.e., the probability of all eigenvalues 
being real for the product XY, where X and Y are from flTSl) . For even N we have 
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while for odd we have 
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Using high-precision numerical computation it was conjectured in [15] that the Meijer 
G-functions in the above expressions lead to values which are rational multiples of tt^. 
We show below that this is indeed the case. 

We set Oi = 02 = 5/2 — j, 6 i = 1 , 62 = ^ + 1 , c = 2 , n = A; — 1 in ffT7|) . which gives 
on little simplification the following remarkable identity 
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This is the central result of this paper which confirms the algebraic and arithmetic 
structure predicted in [15]. Using Eq. (El]), it takes less than a second to evaluate the 
above Meijer G-function with j, k values as large as 1000 in Mathematica m- We have, 
for j, A; = 1,...,5, 
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Interestingly, Mathematica also returns the sum in fET]) in terms of a generalized 
hypergeometric function, which gives 
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Table 1. Exact values and numerical values (6 significant digits for > 1) for 
probability of all eigenvalues real for product of two N x A^-dimensional real 
Gaussian matrices are listed in the second and third columns, respectively. The fourth 
column displays numerical values of the the ratio (4/7r)(p;^^^ ,n-iPn'+i,n+i)/{Pn'^nY 
and supports the leading large N form (7r/4)'^^/^. 


N 

XY 

Pn,n 

4 Pn-i,n-iPn+i,n+i 

Exact 

Numerical value 


1 

1 

1 

— 

2 

71 

7.85398 X 10"^ 

1.01321 

3 

hvr 

4.90874 X 10-1 

1.00500 

4 

20l7r2 

213 

2.42162 X 10-1 

1.00446 

5 

100137r2 

220 

9.42462 X 10-2 

1.00257 

6 

640115857r3 

236 

2.88821 X 10-2 

1.00229 

7 

31625532537A 

247 

6.96751 X 10-® 

1.00156 

8 

8012440011007425A 

269 

1.32219 X 10-® 

1.00142 

9 

3918664131501112628177^ 

284 

1.97341 X 10-^ 

1.00106 

10 

62866533933446109812619543457r5 

2116 

2.31574 X 10-® 

1.00098 

11 

30407079048718892174159408210872571® 

2.13636 X 10-® 

1.00077 

2135 


This result, when combined with leads to yet another interesting identity 

3^2(1 ,3 + k- 1/2, j + k- 1/2; k + G 2j + 2fc - 1; 1) = ~ 

_i^ + _ 

v\3 + k- 1/2)^r(/i + i)r(/i + 2j + k-iy ^ ’ 

We note here that 3 F 2 (ai, 02 , 03 ; A,/^ 2 ; = 3 ^ 2 (ai, 02 , 03 ; A, A; A/(r(A)r(A))- 

With the result fl2lll at our disposal, we can also calculate exact results for 
using (11911 and (|20|) as rational multiples of powers of vr. Exact expressions for this 
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Figure 1. Plots of and the leading large-A^ asymptotic for N up to 

100. The y-axis scale is logarithmic. 


probability for N as large as 100 is obtained in about 5 seconds using Matliematica. 
We list exact results up to = 11 in Table 1, along with the numerical values with six 
significant digits. Also displayed in the third column are the numerical values of the 
ratio (4/7r)(p;^^;^ ,n-iPn+i,n+i)/ ■> which corroborates the large N leading result 
(7r/4)^ [15]. In Fig. [T] we show the plots of exact and the large N leading 

contribution The curves are in good agreement. 

4. Summary and Outlook 

In this work we explored some identities involving Meijer G and hypergeometric 
functions. Using these we proved a recent conjecture of Forrester regarding the algebraic 
and arithmetic structure of Meijer G-function appearing in the study of probability of 
all eigenvalues real for product of two real Gaussian random matrices. Whether similar 
structures exist in the case of product of more than two matrices, remains an open 
question. 
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